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本 


摘 要 : 针对 花 采 授粉 算法 收 你 速度 慢 ， 寻 优 精 度 低 的 缺陷 ， 提 出 基于 折射 原理 的 混合 型 花 采 授粉 算法 (refrHFPA)。 
算法 首先 利用 和 声 搜 索 算法 提升 算法 收敛 速度 ， 然 后 利用 折射 原理 提高 种 群 的 多 样 性 ， 帮 助 算法 跳出 局 部 最 优 ， 提 升 
寻 优 精度 。 实 验 利用 8 个 测试 函数 ， 对 比 其 他 群 智 能 算法 ， 结 果 表 明 refrHFPA 算法 在 收敛 速度 和 寻 优 精度 方面 均 有 显 
著 的 提高 。 
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Hybrid flower pollination algorithm based on refraction principle 
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Abstract: This paper proposed a hybrid flower pollination algorithm (refrHFPA) based on refraction principle for the slow 
convergence rate of flower pollination algorithm and low optimization accuracy. The algorithm firstly used the harmony search 
algorithm to improve the convergence speed of the algorithm, then it used the refraction principle to improve the diversity of 
the population, and helped the algorithm to jump out of the local optimal and improved the accuracy of optimization. It used 
eight test functions to compare other intelligent algorithms, and the results show that the refrHFPA algorithm has a significant 
improvement in convergence Speed and optimization accuracy. 
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0 引言 文献 对 FPA 算法 的 寻 优 能 力 有 一 定 程度 上 的 改进 , 但 仍 存在 寻 
优 精度 低 、 后 期 收敛 速度 慢 等 问题 。 

花 林 授粉 算法 (flower pollination algorithm，FPA) I 是 一 针对 花 条 授粉 算法 速度 和 精度 的 问题 ， 本 文 将 折射 原理 与 

种 新 型 元 启发 式 智能 算法 ， 由 英国 学 者 Yang 于 2012 年 首次 提 ” ”和 声 搜索 算法 引入 到 花 森 授粉 算法 中 ， 提 出 基于 折射 原理 的 泥 

出 。 由 于 FPA 算法 具有 参数 较 少 、 实 现 简 单 、 易 与 其 他 算法 ; 合 型 花茶 授粉 算法 (refrHFPA )。 该 算法 将 和 声 搜索 算法 引入 到 


行 组 合 优化 ， 且 转换 概率 p 可 以 较 好 地 调节 全 局 搜索 与 局 部 搜 。 FPA 算法 中 ， 以 提高 原 算法 的 寻 优 速度 ， 利 用 折射 原理 提高 原 
索 之 间 平 衡 等 优点 ,目前 已 在 函数 优化 领域 得 到 了 应 用 ,例如 ， 种 群 多 样 性 ,使 FPA 算法 跳出 局 部 最 优 ， 以 提高 算法 的 收敛 精 
中 国学 者 Zhou 等 人 外 应 用 改进 的 花 杂 授粉 算法 优化 无 人 海底 ” 度 。 

车 辆 的 路 径 规划 ; 贺 圣 彦 等 人 口 将 该 算法 应 用 到 PID 参数 优化 ; 1 ”基础 理论 

Abdel-Baset 等 人 外 利用 该 算法 进行 比例 优化 。 但 是 花 东 授粉 算 

法 也 具有 仿生 算法 的 典型 缺点 ， 如 易 陷入 局 部 极 值 、 后 期 收 僵 。 1.1 花灯 授粉 算法 

速度 慢 等 问题 为 了 解决 这 些 问 题 , 文献 [5] 将 FPA 与 萤火虫 算 花 末 授粉 算法 是 一 种 模拟 开花 植物 授粉 进程 的 智能 算法 。 
法 相 结 合 ， 在 局 部 搜索 部 分 引入 自 适应 因子 并 自 适应 调整 转换 。 ”授粉 方式 分 为 非 生物 自 花 授 粉 和 生物 异 花 授 粉 。 非 生物 自 花 授 
概率 p。 文 献 [6] 利 用 高 斯 变异 对 全 局 搜索 进行 扰动 增强 种 群 的 ” 粉 是 指 花粉 的 传播 不 需要 生物 作为 传播 者 ， 而 是 通过 风 进 行 伟 
多 样 性 ， 将 Powell 法 引入 局 部 搜索 ， 以 提升 其 局 部 搜索 能 力 。 ” 播 ， 这 种 传播 方式 范围 较 小 ， 被 视 为 局 部 搜索 。 生 物 异 花 授粉 
文献 [7] 将 差分 进化 的 变异 、 交 叉 和 选择 操作 应 用 到 花 条 授粉 算 。 是 两 种 开花 植物 进行 授粉 ， 授 粉 过 程 是 由 蜜蜂 等 传粉 者 进行 传 
法 帮助 算法 跳出 局 部 最 优 ， 以 增强 全 局 搜索 能 力 。FPA 算法 因  ” 播 ， 该 过 程 因 可 在 大 范围 内 进行 传播 ， 被 视 为 全 局 搜索 。 全 局 
提出 时 间 短 ， 迫 切 需要 理论 研究 及 应 用 领域 的 扩展 。 虽 然 以 上 ” 搜索 与 局 部 搜索 之 间 的 转换 通过 转换 概率 p 来 进行 调节 。 现 实 
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生活 中 的 每 一 棵 显 花 植物 都 可 以 产生 数 百 万 甚至 更 多 的 花粉 配 
子 。 为 了 简化 问题 ， 通 常 假 设 每 棵 显 花 植物 只 开 一 水 伦 ， 而 且 


ChinaXiv 合 作 其 
基于 折射 原理 的 混 上 人 


崔 丽 群 ， 等 : 合 型 花 采 授粉 
0) 生 成 一 个 新 的 和 声 。 新 和 声 的 每 一 个 音调 通过 如 下 三 种 


机 理 产 生 : 


每 休 花 只 产生 一 个 花粉 配子 ， 这 代表 每 个 花粉 配子 只 对 应 一 个 


根据 开花 植物 的 授粉 过 程 ， 花 条 授粉 算法 
素 [1]， 

a) 生 物 异 花 授粉 是 携带 花粉 的 传粉 者 进行 全 局 搜索 ， 方 式 
符合 莱 维 飞行 ; 

b) 非 生物 自 花 授粉 被 视 为 局 部 搜索 过 程 ; 

0) 花 的 恒 常 性 被 认为 是 繁殖 概率 ， 这 种 概率 与 所 涉及 的 两 
杂 花 的 相似 性 成 正比 ; 


有 如 下 四 条 规 


一 个 


(a) 学 习 和 声 记 忆 库 。 
Cb ) 若 新 解 第 一 个 变量 以 HMCR 的 概率 选 自 和 声 记 忆 库 ， 
则 按照 公式 (2) 对 音调 进行 微调 ， 反 之 则 在 变量 范围 内 (和 声 
记忆 库 外 ) 随机 选择 一 个 值 作为 新 解 。 
， [x 土 nxBW _ 若 r, < PAR 
| 其 他 名 
其 中 : i、 为 符合 [0,1] 均 匀 分 布 的 随机 数 。 
(c) 随 机 选择 音调 。 
(d) 将 生成 的 新 和 声 和 原 有 和 声 记忆 库 中 的 最 差 值 进行 比 


d) 转换 概率 p e[0,1] ， 受 物理 距离 和 
整个 授粉 活动 中 ，p 值 的 选取 非常 重要 。 
1.2 折射 原理 

光 的 折射 是 指 光 从 一 种 介质 斜 射 入 另 一 种 介质 ， 它 是 自然 
界 中 存在 的 一 种 普遍 的 自然 现象 Bl， 如 图 1 所 示 。 其 中 ol 为 入 
射 朋 ，o 为 折射 角 ，O 为 交界 点 。 光 的 传播 方向 在 不 同 介质 的 
交界 处 发 生 偏 折 。 


其 他 因素 的 影响 ， 在 


入 射 光 反射 光 


空气 


其 他 介质 


图 1 光 的 折射 原理 


斯 涅 耳 通 过 实验 确立 了 折射 定律 ， 即 介质 的 绝对 折射 率 n 
为 入 射 角 和 折射 角 的 正弦 比 外 ， 公 式 为 


n=sinQ /sinQ, (1) 


2 ”改进 的 花灯 授粉 算法 


2.1 基于 和 声 搜索 算法 的 花灯 授粉 算法 
和 声 搜索 算法 (Harmony search, HS) 是 
等 人 09 于 2001 年 提出 的 启发 式 优化 算法 。 由 文献 [10] 可 知 , HS 
算法 由 于 随机 选择 音调 的 作用 使 其 全 局 搜索 能 力 较 强 。 因 此 将 
HS 算法 引入 FPA 算法 中 ， 得 到 基于 和 声 搜 索 算法 的 花灯 授粉 


韩国 学 者 Geem 


PP 


算法 (HSFPA )。HSFPA 算法 通过 将 HS 算法 的 最 优 解 作为 初始 
解 来 提高 算法 的 寻 优 精度 。 和 声 搜索 算法 具体 步骤 如 下 : 


a) 对 HMS (初始 解 个 数 )、HMCR (和 声 记忆 库 考虑 概率 )、 
PAR (和 声 微调 概率 )、 BW (音调 微调 带宽 ) 等 参数 进行 设置 。 

b) 算 法 进行 初始 化 , 随机 生成 HMS 个 和 声 , 形成 和 声 记忆 
库 。 和 声 记忆 库 为 : 


1 L 1 


1 
A oD | fx) 
2 2 
HM=| ” 了 
HMS HMS HMS HMS 
Xl 72 Jo f(x ) 


较 ， 若 优 于 最 差 值 则 对 和 声 记忆 库 进 行 更 新 ， 反 之 则 不 变 。 

(Ce) 判断 是 否 符合 迁 代 结束 条 件 ,符合 则 结束 迁 代 ,不符 
合 则 重复 (c) 和 (d) 直到 符合 条 件 。 
2.2 基于 折射 原理 的 花朵 授粉 算法 
本 文 将 折射 原理 引入 FPA 算法， 以 增加 种 群 的 多 样 性 ， 
帮助 算法 跳出 局 部 最 优 , 提高 算法 的 收敛 速 度 和 寻 优 精度 。 
具体 过 程 如 图 2 所 示 。 将 空气 与 其 他 介质 的 临界 处 看 为 X 轴 ， 
X 是 搜索 区 域 的 上 下 界 ， 取 值 范围 为 a8，b]， 点 Xtes 是 当前 大 
代 的 最 优 解 ，X 为 Xbest 正 上 方 的 入 射 光 点 ，Y 为 点 Xi 的 折射 
点 ，Y1 则 为 点 在 交界 处 的 投影 即 点 Xues 的 折射 点 。 将 OX 
的 距离 设 为 H, OY 的 距离 设 为 Hi, 交界 点 0 为 中 心 点 (atb) 
/2 


斌 


法 线 Y ”折射 光 


图 2 折射 原理 优化 FPA 算法 原理 示意 


由 图 2 可 知 : 
sin@ =((a+b)2- X,)/H (3) 
sin os =(Y -Ca+bD/2/H, (4) 
由 式 (3) (4) 可 知 ,折射 率 n 可 由 式 〈5) 表示 。 
_((a+b)/2-X,.,)/H 


n= sinoQ, /sina, 
! ~ (¥-(a+b)/2)/H, ©) 


设 HHi=f， 可 将 式 〈5) 改 为 
Y=(a+b)/2+((a+b)/2— Xs)/ (fn) (6) 
算法 每 进行 一 次 迭代 ，Xbest 都 通过 式 (6) 进行 折射 ， 得 到 
折射 解 YI， 通 过 判断 Xbest 与 Yi 的 适应 度 值 大 小 ， 对 当前 的 全 
局 最 优 解 进行 更 新 。 
2.3 ”基于 折射 原理 的 混合 型 花朵 授粉 算法 
将 HS 算法 的 最 优 解 作为 FPA 算法 的 初始 解 ， 可 在 一 定 程 
度 上 提升 算法 的 收敛 精度 ， 但 是 HS 算法 的 随机 性 较 大 ， 有 时 
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会 使 算法 陷入 局 部 最 优 。 因 此 利用 折射 原理 增强 种 群 的 多 样 性 ， 
以 利于 算法 跳出 局 部 最 优 ， 提 高 算法 的 寻 优 精度 和 收敛 速度 
改进 后 的 算法 流程 如 图 3 所 示 。 有 具体 实现 步 又 如 下 : 

a) 初 始 化 算法 的 各 参数 。 包 括 p，HMCR,PAR 等 参数 。 

b) 求 出 每 个 解 的 适应 度 值 。 

oc) 利用 HS 算法 求 出 当前 最 优 解 作为 初始 值 。 

d) 当 p>rand (rand 为 [0,1] 均 匀 分 布 的 随机 数 ) 时， 按照 式 
(7) 对 解 进行 更 新 ， 按 照 式 (8) 对 步 长 进行 更 新 ， 反 之 则 按 
照 式 (9) 对 解 进行 更 新 ， 最 后 进行 越界 处 理 。 

X" =X +L(X, — g,) (7) 

其 中 : Xitl、Xi 分 别 是 指向 量 久 ;迭代 t+1 次 和 t 次 后 的 解 ， 8， 
代表 当前 最 优 解 。 


_ A Wsin(TH) 1 


-二 


mi’(s >> so >0) (8) 
A KY 
其 中 : 工 代 表 步 长 ; (4) 是 标准 的 伽 玛 函 数 ，4=1.5。 
和 ”= 和 +E(CX —X,) (9) 
其 中 : Xi 和 Xu 代表 来 自 相同 植 物 不 同 花 的 花粉 ; s 是 [0,1] 均 


匀 分 布 的 随机 数 。 

6) 求 得 当前 最 优 解 和 最 优 解 的 适应 度 值 。 

作对 当前 最 优 解 按 照 式 (6) 进行 处 理 ， 得 到 折射 点 Y1， 然 
点 的 函数 值 进 行 判断 。 若 优 于 当前 最 优 解 ， 则 更 新 全 局 

反之 不 变 。 

的 最 大 迭代 次 数 或 固定 收敛 精 度 ， 

否则 转 至 ¢)。 


后 对 该 
最 优 解 和 最 优 解 的 适应 度 值 ， 

多 判断 是 否 满足 实验 设 
若 满足 则 输出 结果 结束 算法 ， 


3 ”仿真 实验 及 结果 分 析 


为 了 检验 refrHFPA 算法 的 有 效 性 ， 利 用 八 个 测试 函数 在 
MATLAB 的 环境 下 进行 测试 , 并 将 其 测试 结果 与 FPA、HSFPA、 


refrFPA 算法 进行 比较 。 其 中 ， 八 个 测试 函数 如 下 : 
(GD AD = ,he[-5.12,5.12] .i=1,2,3..n, 该 函数 的 理 
论 最 优 值 为 0。 


C2) f(x)= Ys +0.5) ,el-100,100], 1,2,3...n, 
i=l 


该 函数 的 理论 最 优 值 为 0。 


(3) 广 XE[-10,10] ,该 函数 的 


ea 


里 论 最 优 值 为 0。 


(4) fs(7)= 


pr +1 Le x e[—600,600] ， 
二 1,2,3...n, 该 函数 的 理论 最 优 值 为 0。 


， XE[-10,10] ,i=1,2,3..,n, 该 函数 的 


(5) AGO= 祥 js-: 


I 


旱 论 最 优 值 为 0。 


n 2 nn 
(6) f.(x)=20+e—20exp(-0.2 上 ) exp( 工 ycos(2zrx)， 
7 i=1 7 i=1 


a / 八 \ // 
ee 山 naxIV 合 | | FE 
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XE[-32,32] ,这 1.2,3.n, 该 函数 的 理论 最 优 值 为 0。 


的 理论 最 优 值 为 0。 


x E[-10,10] ,二 1.2,3.n, 该 函数 


(8) 大 CD= >do7 se[100100],i-L23.m 该 


数 的 理论 最 优 值 为 0。 


计算 每 个 解 的 适应 度 值 


HS 算法 求 出 当前 
最 优 解 作为 初始 值 


按照 式 (7) 对 解 进行 更 新 


按照 式 (8) 对 步 长 进行 更 新 


新 全 局 最 优 解 和 最 优 
的 适应 度 值 


按照 式 (9) 对 解 进行 更 


获取 折射 点 Yi 


上 点 优 于 当 
前 最 优 解 ? 


i 全 局 最 优 解 和 最 优 
0 


图 3 ”refrHFPA 算法 流程 


3.1 参数 f 的 分 析 与 选取 

参数 f 的 取 值 直接 影响 折射 解 在 X 轴 上 的 位 置 ， 取 了 在 预 
估 范 围 [0,1] 内 , 对 单 峰 函数 访 和 万、 多 峰 函 数 f 和 无 独立 运行 
20 次 。 判 断 f 取 何 值 时 可 以 有 效 帮 助 算 法 跳出 局 部 最 优 。 本 实 
验 设置 种 群 数 为 20， 最 大 友 代 次 数 为 300，refrHFPA 算法 参数 
中 p=0.201, 和 =1.5, HMCR=0.95, PAR=0.3, BW=0.01, n=28]。 
实验 结果 如 表 1 所 示 。 
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表 1 f 取 不 同 值 时 refrHFPA 算法 的 性 能 分 析 
函数 维 数 f 最 优 值 平均 值 最 差 值 
0.2 2.492e-01 $5.4612e-01 9.545e-01 
0.4 1.675e-01 7.1101e-01 1.2157 
f 30 0.6 3.Se-03 2.3523e-02 1.31Se-01 
0.8 6.1e-04 3.7303e-03 8.6e-03 
1 6.62e-04 1.6072e-03 3.1e-03 
0.2 1.6512e+02 2.7669e+02 3.6554e+02 
0.4 1.0628e+02 3.0286e+02 7.2461e+02 
万 30 0.6 8.3564 1.5407e+01 2.4S4Se+01 
0.8 3.0175 8.2552 4.5821e+01 
1 2.8153 4.9436 7.2475 
O02 2.4283 3.82006 5.0583 
0.4 1.7756 3.71579 6.7429 
fa 30 0.6 1.0333 1.06872 1.1426 
0.8 3.353e-01 8.3804e-01 1.0272 
1 2.949e-01 7.7102e-01 1.0487 
0.2 5.2427 6.94923 1.0569e+01 
0.4 4.4971 7.02611 1.1209e+01 
fo 30 0.6 8.233e-01 2.49892 4.2241 
0.8 3.862e-01 9.3248e-01 1.6789 
1 1.898e-01 S.7248e-01 1.253 


优 结果 的 影响 ， 实 验 采用 与 ， 收 敛 情况 如 图 


fl 函 数 收敛 曲线 


步 确定 参数 


表 1 可 以 看 出 ， 参 数 f 取 值 [0.6,1] 相 比 [0, 0.4] 寻 优 结果 
有 1~3 个 数量 级 的 提高 。 为 了 进 


取 值 对 算法 寻 


f4 函 数 


党 


NGY 


4 所 示 。 


by 


| 


图 4 不 同 f 取 值 测 试 函 数 的 收敛 曲线 


照 式 〈10) 进行 取 值 ， 初 始 值 设 为 0.8。 


3.2 测试 方案 
本 文 共 设计 五 利 
迭代 次 数 ， 通 过 八 利 


实验 方案 : 第 


FE 


refrHFPA 算法 在 低 维 天 


卉 下 的 寻 


f =0.8+0.2xrand() 


4 可 见 , f 值 在 [0.8,1] 内 时 收敛 速度 更 快 ， 因 此 将 f 按 


(10) 


种 是 固定 种 群 数量 及 
上 试 函 数 测试 FPA、HSFPA、refrFPA 及 


入 癌 - 


维 的 实验 环境 , 判断 refrHFPA 算法 在 高 维 


人 能 力 ; 第 二 种 是 设 定 局 


环境 下 是 否 可 以 得 到 


较 好 的 收敛 精度 ; 


第 三 种 是 


回 定 | 


妇 仇 精度， 对 比 这 


四 种 算法 的 


最 小 收敛 次 数 、 平 均 收 敛 次 数 和 收敛 率 ; 第 四 种 是 refrHFPA 算 


法 在 不 同 维度 下 ， 风 


是 将 refrHFPA 算法 与 近 两 各 


本 文 设 


证 本 文 算法 的 可 行 性 
的 实验 参数 如 下 : 


tT 


| 断 其 是 否 具有 较 低 的 时 间 复 杂 
改进 后 的 FPA 算法 进行 对 比 ， 验 


度 ; 第 五 种 


FPA 算法 


BW=0.1。 


refrFPA 算法 中 p=0.2[ 


1 
| 


上 岩 丽 群 ， 等 : 基于 匠 射 原理 的 混 


转换 概率 p=0.801， 和 =1.5。 
HSFPA 算法 中 p=0.8 


，A=1.3，HMCR=0.9，PAR=0.3， 


1, MA=1.5,，n=2[8],， 人 0.8, 。 


生 能 分 析 


迭 


refrHFPA 算法 参数 设置 与 小 节 3.1 相同 。 
3.2.1 国定 选 代 次 数 ， 
1) 低 维 下 的 性 能 分 析 
本 文 实验 方案 设置 种 群 数 为 20， 


民 次 数 为 1000， 态 维 


数 为 10, 其 余 函 数 维 数 为 30。 每 种 算法 对 每 个 测试 函数 均 独 立 
运行 50 次 ， 计 算 其 最 优 值 、 最 差 值 、 平 均值 和 寻 优 成 功率 。 
其 中 寻 优 成 功 的 标准 是 实际 最 优 值 与 理论 最 优 值 的 差 值 小 于 
0.001。 实 验 结果 如 表 2 所 示 。 
表 2 refrHFPA 算法 在 低 维 上 的 性 能 分 析 
函数 ”算法 平均 值 最 差 值 最 优 值 成 功率 
2.74751 4.483 1.5283 0 
HSFPA 0.06378 0.13458 0.01626 0 
n refrFPA 1.8486e-08 4.4378e-08 3.426e-09 100% 
refrHFPA 9.4645e-12 1.6008e-10 7.514le-15 100% 
EpA 12752e+03 3.3109e+03 1.0687e+02 0 
HSFPA 2.9799e+01 4.4754e+01 1.1784e+01 0 
h refrFPA 5.696e-04 1.1301e-03 2.1231e-04 92.8% 
refrHFPA 2.7042e-05 1.3952e-04 2.5498e-06 100% 
ey 1.48352 2.429 5.389e-01 0 
HSFPA 5.8429e-03 3.06e-02 1.9e-03 0 
h refrFPA 9.2467e-01 2.5223 6.5e-03 0 
refrHFPA 9.1811e-09 8.6994e-08 3.3573e-10 100% 
EpA 14417e+01 1.1886e+02 4.77 0 
HSFPA -1.37221 1.9232 1.1005 0 
h refrFPA 5.0578e-03 4.34e-02 4.6823e-06 50% 
refrHFPA 1.4265e-10 5.89e-09 5.9064e-14 100% 
EpA 131l2e+01 1.8652e+01 9.1526 0 
HSFPA 1.01672 1.4035 0.6329 0 
上 refrFPA 2.3643e-03 3.5e-03 1.8e-03 0 
refrHFPA S$.1908e-04 9.34Sle-04 2.6926e-04 100% 
A 5.81155 7.9429 2.8798 0 
HSFPA 2.95162 4.4069 2.036 0 
refrFPA 7.8168e-02 7.838e-01 8.5053e-04 7% 
refrHFPA 2.4028e-05 4.0876e-05 1.8607e-05 100% 
FPA 1.495le+1l 2.2939etl 9.4265 0 
HSFPA 1.58679 2.6311 4.712e-01 0 
h refrFPA 7.8514e-05 1.6539e-04 2.9294e-05 100% 
refrHFPA 1.248e-07 1.5261le-06 4.3813e-09 100% 
EpA 1.0406e+05 2.8598e+05 3.392e+04 0 
HSFPA 4.7374e+04 7.2367e+04 1.9956e+04 0 
h refrFPA 1.7043e-02 5.08e-02 1.3e-03 0 
refrHFPA 2.6904e-07 2.9884e-06 ”9.7789e-11 100% 


法 的 最 全 
算法 ; 


值 、 最 差 值 及 了 
除 外，HSFPA 较 于 
数量 级 的 提高 ， 除 户外 ， 对 了 


从 表 2 可 以 看 出 ， 针 对 本 文选 取 的 测试 函数 ，refrHFPA 算 
均值 均 优 于 FPA、HSFPA 和 refrFPA 
FPA 算法 ， 寻 优 结果 均 
不同 的 测试 函数 refrFPA 算法 较 


1~2 个 


于 原 算 法 均 有 4~8 个 数量 级 的 提高 。 此 外 ，refrHFPA 算法 对 于 
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录用 入 上 岩 丽 群 ， 等 : 基于 une 和 采 授 粉 算 法 
这 八 种 测试 函数 的 寻 优 率 达 到 100%, 并 且 最 优 值 、 最 差 值 和 平 表 3 refrHFPA 算法 在 高 维 上 的 性 能 分 析 
均值 也 远 高 于 其 他 三 种 算法 ,为 了 更 直观 地 比较 refrHFPA 算法 函数 维 数 ”算法 平均 值 最 差 值 最 优 值 ” 成 功率 
与 其 他 三 种 算法 寻 优 能 力 ， 本 文 对 于 八 个 测试 本数 的 收 全 情况 jo mA S40109 S383et03 28358e103 0 
如 图 5 所 示 。 2 refrHFPA -1.5171e-07 3.586e-07 4.5143e-08 100% 
FPA 5.1105e+01 8.8088e+01 3.1268e+01 0 
100 
f1 函 数 收敛 曲线 f2 函 数 收 敛 曲 线 A 4 refrHFPA 2.5373e-13 3.6189e-12 3.3307e-16 100% 
140 80000 f FPA 5.6286 9.3972 2.8912 0 
20 100 
地 60000 6 refrHFPA 1.8174e-05 1.8273e-05 1.8172e-05 100% 
加 | 0 FPA 4.7518e+07 8.7437e+07 1.9636e+07 0 
2 当 20009 及 refrHFPA 5.6762e-06 1.9884e-04 7.2084e-12 100% 
0 Ee 0 i 
1 201 401 601 801 3 201 401 601 801 ,i a 
0 ee a ee 表 3 表明， 在 高 维 环境 下 refrHFPA 与 FPA 算法 相 比 ， 在 
~- refrFPA 一 一 reffHFFA | refrFPA 一 一 refrHFPA 寻 优 精度 上 有 5~15 个 数量 级 的 提高 。 因此 refrHFPA 算法 可 以 
f3 函 数 收 敛 曲线 全 函数 收敛 曲线 跳出 局 部 最 优 ， 不 易 随 着 维 数 的 增加 而 陷入 “ 维 灾难 ”。 
2.00Ex01 - 2000 3.2.2 定 收 化 精度 性 生 能 分析 
{ 800 
:aoea 加 收敛 率 是 寻 优 成 功 的 次 数 占 实验 总 次 数 的 百分比 。 为 了 验 
PT 400 证 refrHFPA 算法 的 收敛 性 , 在 固定 精度 的 前 提 下 ,对 每 个 测试 
5.00E+00 
Es ~ 发 函数 独立 运行 50 次 , 将 其 最 小 收敛 代数 、 平均 收敛 代数 、 收敛 
401 601 801 a 率 与 FPA、HSFPA、refrFPA 算法 相 比 较 。 实 验 中 四 种 算法 的 种 
—— FPA HSFPA 
ee refrFPA 一 一 refrHFPA 人 群 数 均 为 20， 维 数 设 定 为 10， 最 大 人 迭代 次 数 为 2 500。 实 验 结 
条 函 数 收敛 曲线 果 如 表 4 所 示 。 
5 函数 收敛 曲线 
人 表 4 四 种 算法 在 固定 收敛 精度 下 的 性 能 分 析 


函数 维 数 固定 精度 算法 最 小 收敛 代数 平均 收敛 代数 ”收敛 率 


oT FPA 961 1225.58 100% 

40 

20 3 一 一 一 一 -0 HSFPA 582 685.97 100% 
0 - i 10 1.0e-05 


1 201 401 601 801 


1 201 401 601 801 refrFPA 224 290.43 100% 
一 一 FpPA HSFPA Fn i refrHFPA 161 216.18 100% 
-refrFPA 一 一 refrhFPA | ---- refrFPA 一 一 refrHFPA 

FPA 1760 2059.17 82% 
函数 收敛 曲 ee HSFPA 1076 1182.36 100% 
全 函数 收敛 曲线 f 函 数 收 敛 曲 线 f 10 1.0e-07 
160 2 refrFPA 633 666.14 100% 
| 2.50E:09 refrHFPA SS8 601.08 100% 

2.00E+09 | FPA fai ail 0 
60 | HSFPA 927 1193.71 100% 

40 1.00E*09 fs 10 1.0e-03 ee 

a Tr 00go0s | refrFPA fai ail 0 

a ee 也 refrHFPA 325 397.36 100% 
201 401 601 801 
FPA ep — mA HSFPA ----- refrFPA 一 一 refrHFPA 全 fa fail d 
mae refrFPA 一 一 refrHFPA HSFPA fai fail 0 
六 10 1.0e-07 

图 5 四 种 算法 对 于 测试 函数 的 收敛 曲线 攻 fai ail 
refrHFPA 202 594.38 100% 
从 图 5 可 以 看 出 , HSFPA 算法 在 对 测试 函数 方 -万 进行 测 FPA 1076 1165.2 100% 

未 i a Ey HSFPA 872 959.79 1009 
试 时 陷入 局 部 最 优 ， 但 refrHFPA 算法 却 能 很 好 地 跳出 局 部 最 大 io om ， ， 。 
5 refrFPA 437 455.43 100% 

、 十 全 站 右 施 的 丧 : 半 尾 尘 

优 ， 说 明 引入 折射 定理 可 以 有 效 的 增加 种 群 多 样 性 ， 帮 助 算法 SEE 
跳出 局 部 最 优 ， 提 升 算法 的 收敛 速度 和 寻 优 精度 。 FPA 1542 2056.3 86% 
Si HSFPA 910 1045.07 100? 
2) 高 维 下 的 性 能 分 析 fo 0 0 
yo , 6 refrFPA 500 730.29 100% 
目前 处 理 的 很 多 数据 多 为 高 维 数据 ， 如 多 媒体 图 形 图 像 和 ea 
视频 数据 等 , 本 文 算法 在 高 维 环境 下 进行 实验 。 维 数 设 为 100， FPA 1856 2115.23 62% 
友 代 次 数 设 为 4000， 其 余 参 数 不 变 。 对 于 每 个 测试 函数 ， 将 原 f 10 10e.03 SEEA A 和 So 
多 报 粮 税法 上 木 立 担 多 授 粮 位 法 均 独 立 字 验 gp 7 refrFPA 258 344.36 100% 
花 杂 授粉 算法 与 本 文 提出 的 花 杂 授粉 算法 均 独立 实验 50 次 ， = 
若 实 际 最 优 值 与 理论 最 优 值 的 差 值 小 于 0.001 时 代表 寻 优 成 功 。 FPA 2271 2479.1 20% 
实验 结果 如 表 3 所 示 。 f de HSFPA 1486 1611.21 100% 
8 refrFPA 465 606 100% 


refrHFPA 450 582.6 100% 
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从 表 4 可 以 看 出 ， 对 于 这 八 个 测试 函数 ，HSFPA、reftFPA ”中 的 部 分 函数 。 本 实验 方案 的 参数 设置 与 参考 文献 相同 ， 
与 refrHFPA 算法 无 论 是 最 小 收敛 代数 还 是 平均 收敛 代数 均 优 。 GMPFPA 和 DEFPA 算法 的 实验 数据 均 来 自 于 参考 文献 ， 结 果 
于 原 算法 。 由 此 可 见 HS 算法 与 折射 定理 的 引进 对 于 FPA 算法 如 表 6 和 7 所 示 。 
均 是 有 效 的。 此 外 ，refrHFPA 算法 的 最 小 收敛 代数 和 平均 收敛 表 6 reftHFPA 算法 与 GMPFPA 算法 优化 性 能 对 上 
代数 在 四 种 算法 中 最 小 ， 且 收敛 率 均 为 100%。 因 此 refrHFPA 函数 ” 维 数 。 ”算法 平均 值 
算法 的 收敛 性 要 优 于 其 他 三 种 算法 ， 并 具有 有 效 性 。 


N| 


CL 


t 


最 差 值 最 优 值 


GMPFPA 1.4323e-18 1.5429e-17 1.0003e-30 
3.2.3 refHFPA 算法 的 时 间 复 杂 度 分 析 和 refrHSFPA 7.5053e-30 4.2212e-29 3.4197e-33 
本 文 将 refrHFPA 与 FPA 算法 在 相同 环境 下 ， 通 过 运行 时 GMPFPA 1.1166e-13 1.3006e-12 7.886e-13 
间 来 观察 本 文 算法 是 否 满足 时 间 复杂 度 较 低 的 要 求 。 实 验 中 种 J 
群 数量 设置 为 10， 和 迭代 次 数 为 300， 维 数 分 别 设 为 10、100、 GMPFPA 1.388le-10 1.467e-09 4.4409e-15 
200， 每 种 算法 对 于 每 种 测试 函数 均 独 立 运行 20 次 。refrHFPA j6 30 refrHSFPA 5.7731le-15 1.0658e-14 3.5527e-15 
与 FPA 算法 在 不 同 维度 下 的 最 小 运行 时 间 与 平均 运行 时 间 如 GMPFPA 1.501le-16 2.3621e-15 2.7323e-21 
表 5 所 示 。 J refrHSFPA 6.0552e-19 2.7656e-18 2.2322e-20 
表 5 refrHFPA 算法 时 间 复 杂 度 分 析 
函数 维 数 算法 最 优 值 平均 值 表 7 refrHFPA 算法 与 DEFPA 算法 优化 性 能 对 比 
出 FPA 0.302 0.31624 函数 。 维 数 算法 平均 值 最 差 值 最 优 值 
refrHFPA 0.902 0.94136 DEFPA 2.6347e-17 1.0206e-16 7.4022e-18 
FPA 0.603 0.62993 J refrHSFPA 1.8415e-32 6.587e-31 1.0809e-34 
1 refrHFPA 1.662 1.70079 DEFPA 5.0093e-15 1.4322e-14 9.992e-16 
FPA 0.929 0.966 坟 各 refrHSFPA 0 0 0 
refrHFPA 2.453 2.53729 DEFPA 3.359e-08 6.9112e-08 1.4435e-08 
FPA 0.339 0.3765 j6 2 refrHSFPA 7.993Se-15 2.1316e-14 3.5527e-15 
refrHFPA 0.903 1.30135 DEFPA 9.216le-31 2.2095e-30 3.66e-31 
FPA 0.703 0.74143 J refrHSFPA 4.6329e-19 1.8917e-18 1.5876e-20 
有 加 refrHFPA 1.662 2.7271 
FPA 1081 1 12068 由 表 6 可 以 看 出 ，refrHFPA 算法 相 较 于 GMPFPA 算法 ， 
200 nk 本 二 有 单 峰 函 数 有 12 个 数量 级 的 提高 ， 多 峰 函 数 有 3~5 个 数量 级 的 
¢ FPA 034 0.35907 是 高 。 由 表 7 可 以 看 出 ， 除 户外，refrHFPA 算法 较 于 DEFPA 
一 ee 3 a 算法 , 在 寻 优 结果 上 均 有 6~16 个 数量 级 的 提高 。 因 此 , 在 实验 
SR Dred 0 条 件 相同 的 情况 下 ， 无 论 是 单 峰 函数 还 是 多 峰 函 数 ，refrHFPA 
f ee refrHFPA 2.458 2.60107 算法 寻 优 结果 优 于 参考 文献 [6,7] 算 法 。 
本 FPA 1.079 1.15293 4 ”结束 语 
refrHFPA 3.838 4.02456 
和 i Bee 花 东 授粉 算法 作为 新 型 群 智能 算法 ， 也 与 其 他 智能 算法 一 
.0 fA ye 0 样 存在 收敛 速度 慢 、 寻 优 精度 低 的 缺点 。 针 对 这 些 不 足 ， 本 文 
Ne ee B20 将 HS 算法 引入 花茶 授粉 算法 ， 将 其 最 优 解 作为 算法 欠 代 的 初 
广 By pe 1 5 台 解 ， 以 提高 初始 解 的 质量 。 利 用 折射 定理 对 当前 最 优 解 进行 
人 Ds a 折射 ， 以 提升 算法 的 全 局 搜索 能 力 ， 提 高 寻 优 精度 。 实 验 结果 
eu DA 2 ed 表明 ,本文 算法 与 原 算法 和 其 他 改进 算法 相 比 ,收敛 速度 较 快 、 
最 优 解 质量 较 高 ， 具 有 一 定 的 可 行 性 和 有 效 性 。 以 后 ， 将 在 花 
从 表 5 可 以 看 出 ，refrHFPA 与 FPA 算法 相 比 运行 时 间 略 。 杂 授 粉 算法 的 参数 值 设 置 、 推 广 算法 应 用 领域 等 方面 做 进一步 


长 , 但 refrHFPA 算法 无 论 在 低 维 还 是 高 维 环境 下 的 寻 优 精度 都 。 研究 。 


远 高 于 FPA 算法 ， 而 且 refrHFPA 算法 的 最 小 运行 时 间 与 平均 参考 文献 
fe 2 六 人 : 
运行 时 间 仍 然 较 短 ， 所 以 可 以 看 出 refrHFPA 算法 是 可 行 的 。 四 
3.2.4 改进 智能 算法 比较 实验 [1] Yang Xinshe. Flower pollination algorithm for global optimization [C]// 
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201804.02072V1 


chinaXiv 


录用 稿 


[3] 


5] 


6] 


7] 


8] 


Natural Computation. 2012: 240-249. 

Zhou Yongquan, Wang Rui. An improved flower pollination algorithm for 
optimal unmanned undersea vehicle path planning problem [J]. International 
Journal of Pattern Recognition and Artificial Intelligence, 2016, 30 (4): 
1659010. 

贺 圣 彦 ， 草 中 清 ， 余 胜 威 . 基于 花 授粉 算法 的 PID 参数 优化 [四 . 计算 
机 工程 与 应 用 , 2016, 52 (17): 59-62. 

Abdel-Baset M, Hezam I M. An improved flower pollination algorithm for 
ratios optimization problems [M]// Applied Mathematics and Information 
Sciences Letters. 2015: 83-91. 

卞 京 红 , 贺 兴 时 , 杨 新 社 . 基于 营 火 虫 算法 的 自 适 应 花 授粉 优化 算法 
四 . 计算 机 工程 与 应 用 , 2016, 52 (21): 162-167. 

肖 辉 辉 ， 万 常 选 ， 段 艳 明 ,等 . 融合 高 斯 变异 和 Powell 法 的 花 采 授粉 优 
化 算法 [J]. 计算 机 科学 与 探索 , 2017, 11 (3): 478-490. 
肖 辉 辉 , 万 常 选 , 段 艳 明 . 一 种 改进 的 新 型 元 启发 式 花 采 
计算 机 应 用 研究 , 2016, 33 (1): 126-131. 

邵 胸 ， 吴 志 健 , 周 炫 余 , 等 基于 折射 原理 反 向 学 习 模 型 的 改进 粒子 群 
算法 [J]. 电子 学 报 , 2015, 43 (11): 2137-2144. 


授粉 算法 加]. 


期 刊 


花 杀 授粉 算法 


ChinaXiy 合 1 


瞧 丽 群 ， 等 : 基于 折射 原理 的 混合 型 


[9] Griffiths, David J. Introduction to Electrodynamics [M]. [S. 1. ] : PrentIce 
Hall, 1998. 386-389. 

[10] Geem Z W, Kim J H, Loganathan G V. A new heuristic optimization 
algorithm: harmony search [J]. Simulation, 2001, 76 (2): 60-68. 

[11] Draa A. On the performances of the flower pollination algorithm-qualitative 
and quantitative analyses [J]. Elsevier Science Publishers B. V. , 2015, 34 
(C): 349-371. 

[12] 曲 良 东 , 何 登 旭 ,， 黄 勇 . 自 适应 改进 和 声 一 一 单纯 形 进 化 算法 研究 [J]. 
计算 机 应 用 研究 , 2013, 30 (3): 676-678. 

[13] Zhao S$ Z, Suganthan P N, Pan Q K, et al. Dynamic multi-swarm particle 
swarm optimizer with harmony search [J]. Expert Systems with Applications, 
2011, 38 (4): 3735-3742. 

[14] He X, Yang Xinshe, Karamanoglu M, et al. Global convergence analysis of 
the flower pollination algorithm: a discrete-time markov chain approach [J]. 
Procedia Computer Science, 2017, 108: 1354-1363. 

[15] Zaem D, Mezioud C, Draa A. On the efficiency of the binary flower 
pollination algorithm: application on the antenna positioning problem [j]. 


Applied Soft Computing, 2016, 47 (C): 395-414. 


